UCI-TR-2008-16, MIFP-08-09, NSF-KITP-08-63, arXiv:0805.1229 



Tachyon Mediated Non-Gaussianity 



OO 
O 

o 

(N 

1—5 

OS 
(N 



Or 

cd : 

> 
o> 

(N 

O 

OO 

O 



X 



Bhaskar DuttaEl and Louis LebloncQ 

Department of Physics, Texas A&M University, College Station, TX 77843, USA 

Jason Kumai@ 

Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA 

We describe a general scenario where primordial non-Gaussian curvature perturbations are gener- 
ated in models with extra scalar fields. The extra scalars communicate to the inflaton sector mainly 
through the tachyonic (waterfall) field condensing at the end of hybrid inflation. These models 
can yield significant non-Gaussianity of the local shape, and both signs of the bispectrum can be 
obtained. These models have cosmic strings and a nearly flat power spectrum, which together have 
been recently shown to be a good fit to WMAP data. We illustrate with a model of inflation inspired 
from intersecting brane models. 
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Introduction. Single field slow-roll inflation generically 
predicts a near flat Gaussian spectrum of fluctuations 
in the curvature of spacetime. Gaussianity is a conse- 
quence of the slow-rolling conditions which require the 
inflaton to behave like a free field with negligible self- 
interactions. We can parameterize the non-Gaussianity 
(NG) in the curvature perturbation (denoted £ here) us- 
ing the so called Jml formalism: 



C(X, t) = (Gauss + g/iVlXCc an -S'.s 



(Caauss)) ) (1) 



where (causs is the Gaussian piece of the curvature per- 
turbation. This ansatz for the curvature perturbations 
leads to the so-called local shape for the bispectrum in 
Fourier space. 

The recent five year analysis from Wilkinson Mi- 
crowave Anisotropy Probe (WMAP5) [l| shows —9 < 
f NL < 111 while Yadav and Wandelt (using WMAP3 
data) have claimed a detection of primordial NG at above 
99.5% confidence level 0. It is expected that the upcom- 
ing Planck experiment should be able to detect /jyi, of 
order unity. In the light of all these new and future mea- 
surements, it is very important to investigate the origin of 
large non-Gaussianity in models of inflation. It is also im- 
portant to search for predictions that are associated with 
a larger Jnl for other future observations, e.g., the shape 
of the bispectrum, running of /ml, etc. The best bound 
on Jnl given above is for the local shape. Another well 
known shape, the equilateral shape, has also been con- 
strained (using WMAP3 data) to be -256 < fjf < 332 
Q while many other shapes have not yet been searched 
for individually. 

For single field smooth slow-roll models with canonical 
kinetic terms, //vx is of order of the slow-roll parame- 
ters and is less than 1 In this paper, we investigate a 
natural way to induce large non-Gaussianity from a ubiq- 
uitous feature of string theory and models of physics be- 
yond the standard model: extra scalars. In the following, 
we will refer to the fluctuation along the inflationary tra- 



jectory the adiabatic mode while orthogonal fluctuations 
are entropy modes. The entropy modes fluctuations are 
generically non-Gaussian since in general their potential 
is not flat. Given a mechanism to feed the entropy modes 
fluctuations into curvature perturbations, the curvature 
spectrum will also be non-Gaussian. In addition to this 
intrinsic NG of the entropy modes, the non-linearity of 
the transfer mechanism between entropy modes and cur- 
vature perturbations will also generate significant NG. 

The scenario we present is general but can be easily re- 
alized by intersecting brane models (IBMs) in string the- 
ory The IBMs are well-motivated since in these sce- 
narios not only the standard model can be constructed, 
but also the fermion mass hierarchies can be understood. 
From the effective field theory point of view, this scenario 
gives a multi-fields generalization of D-term inflation [|| . 
We will show that non-Gaussianity can be generated at 
the end of inflation in this type of models from two main 
different sources. The first source of NG is an "intrinsic" 
component coming from self interactions of the entropy 
modes. The second source is coming from the non-linear 
terms that arises when one transfers the entropy pertur- 
bations into curvature perturbations. 

The intrinsic component of non-Gaussianity in our 
work is similar to the ref. Q (building on previous work 
Q) which also considers D-term inflation. The non-linear 
local piece on the other hand has been studied previously 
for multiple field hybrid inflation in Q or in the context of 
modulated reheating [10| . But we find that in more gen- 
eral models, both contributions can appear and can be 
significant. Moreover, we argue that loop-contributions 
to the bispectrum can also be important. Thus, even 
simple, well-motivated models can yield a much more 
complicated picture than previously discussed. 

The intrinsic and non-linear non-Gaussian contribu- 
tions both have local shape with opposite signs that de- 
pend on the hidden-sector charge of the extra fields and 
on whether the dynamics of the entropy mode is domi- 
nated by the gauge sector or the Yukawa sector of the 
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theory. Thus, we find in this setup that the sign of Jnl 
is not a stochastic variable but rather fixed by the under- 
lying physics of the system. We find that the best regime 
to obtain large NG (10 ~ 100) is the regime of D-term 
inflation where we produce cosmic strings with tension 
G/i ~ 10~ 7 and where n s ~ 1. This is a good fit to the 
WMAP data as recently shown in 11]. 

The 2-pt Correlation Function. In this work, we will 
be interested in a model of hybrid inflation containing the 
usual inflaton <f> and the tachyonic (waterfall) field T but 
extended with an extra scalar field x- During inflation, 
the tachyon has a positive mass and its VEV is pinned 
at zero. Below some critical value (<p c ) of the inflaton, 
the mass of the tachyon becomes negative and it con- 
denses, ending inflation. We will assume that the x ne ld 
has low mass, fluctuates during inflation but does not di- 
rectly couple to the inflaton field. Instead its only effect 
is to modulate the mass of the tachyon and to change the 
critical value at which inflation ends. In some ways, the 
tachyon acts as a messenger between the hidden sector 
field x and the inflaton field 0. The action for this system 
is of the following form 

s = \ J VdlMpR W) 2 - (dTf - {d X f - 2V] , 
V = V iai (<t>) + V hid ( X ) + V mess (<t>, X , T)] . (2) 

During inflation, the tachyonic field has vanishing VEV 
and the messenger part of the potential V mC ss vanishes 
as well (all coupling in this part of the potential involves 
the field T). Assuming that the hidden sector potential 
contribute a negligible amount to the vaccuum energy, 
the Hubble scale during inflation is equal to H 2 = . 

This model is of the form considered in @, 0] • Using 
the SN formalism [l2| . we can compute the curvature 
perturbation C(x,t) = SN where N is the number of 
cfolds. Since the inflationary potential is only a function 
of 4>, the field x only affects the critical value at which 
inflation ends. The fluctuation in the number of efolds 
when <ft — * 4> + &4> and X ~~ * X + $X is 



SN 



H 



IHd 2 



Hd4> Cj . 
~T~^~"X 

<t> ox 



f - (sx 2 )) 



(3) 



2 d X 



where the subscript * means evaluated at horizon cross- 
ing. We have kept only the leading term at second order 
and we substracted a number to ensure that (SN) = 0. 
We will find that non-Gaussianities are generated both 
at linear order in the variable Sx as well as from the 
quadratic term. In curvaton models [H (or more re- 
cently in ekpyrotic models [14]), NG is generated from 
these second order terms. The validity of the SN expan- 
sion to higher orders and the NG that are generated from 



these terms were first discussed in 15j. We define 



d(p c 
dx 



(4) 



and this function encodes the coupling between x and 4> 
through the tachyonic field. To evaluate the 2-pt function 
we use the standard result that (ScfiStfi) ~ H 2 /An 2 . Note 
that both S<f> and Sx are not exactly constant outside of 
the horizon and they will decay exponentially because of 
their mass terms with S\\f ~ where 



- /*/ v x Hdt 



-V X Ne 



(5) 



and we used the following definitions for the slow-roll 
parameters. 
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■AH 2 



(6) 



Since we need to evaluate Sx at the end of inflation, sig- 
nificant non-Gaussianity will only arise if this damping is 
not too large, requiring r\ x to be small. However, V iX need 
not be particularly small, which distinguishes x from the 
inflaton. The 2-pt curvature correlation function is there- 
fore 



c(*i)c(fe: 



(2<K) 6 5 6 (k l + k 2 )2<K 2 V^k- 
H 2 ( 1 



8tt 2 A/ 



J 2 K 2 



(7) 



where / means evaluated at the end of inflation. In the 
last expression, we have not included loop contributions 
coming from the higher order terms in the SN. They 
can be made small as we will discuss in more detail later 
on. The spectral index is easily computed from n s — 1 = 

"TEX ' ^ ne f° rmu la simplifies in the limit where e* « e/, 
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1 + 7 2 K' 



(-66* + 2r^ + 7 2 K 2 (-2e, + 2r) x )) 



3-pt Correlation Function. The 3-pt correlation func- 
tion can also be calculated using the same formalism al- 
though in this case we can neglect the part coming from 
the inflaton cj) since it is going to be suppressed by power 
of the slow-roll parameters [4j . The field x on the other 
hand can be non-Gaussian for a general potential Vhid(x)i 
and this non-Gaussianity can be transferred to the curva- 
ture spectrum. We will denote by (C 3 ) int the piece of the 
bispectrum which is intrinsic to the field Sx and arises 
directly from its higher-order self-interactions. This non- 
Gaussianity in Sx is transferred to the SN spectrum by 
the linear terms in eq. ([3]), yielding 



H 



« 3 > int = T T <^ 3 > 



(8) 
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One can calculate this term using the in-in formalism [4| 
with an interaction Hamiltonian simply given by Hi = 
f dPx-iJ^gVhidM- This computation was first done in 
111] (see also [IT 
result is: 



for more recent derivations). The 



£fei ^k 2 Cfc 3 
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(27r) 3 ,5 3 (^) 



H 2 V, 



N, 



8V2M 3 pl e}Ukf 



i<j 



-h 



(9) 



where j e is the Eulcr number and we have defined 
k t = J2^i- The number of efolds is N e = — log( a ) 
and we have used the fact that in the slow-roll regime 



Since the number of efolds from horizon 



h _ _ 

4, \/2eM p 

crossing is around 60 from observations, we see that the 
dominant term is I3 ~ ^ ^ kf ■ This is exactly the lo- 
cal shape obtained from Eq. fTJ) (neglecting the weak 
momentum dependence in N e ). Indeed, Fourier trans- 
forming ( in {l}, one can compute the local 3-pt function 
to be 



(10) 



Comparing the last two expressions we find the following 
intrinsic contribution to the $nl- 



f'mt 
JNL 



5 iV e Af p 7 3 K 6 



9^2 



IPe 3 / 2 



1 



V, 



.(11) 



where we use the subscript "int" to remind the reader 
that this NG is coming from self-interactions of the en- 
tropy modes. This kind of intrinsic NG occurs in models 
of inflation with non-trivial kinetic terms and they usu- 
ally have very distinctive shapes in momentum space (see 
for example [H, Qjl]). In our case, the 3-pt function is 
dominated by the term proportional to the number of 
efolds which originates from the classical evolution of the 
field outside of the horizon and it therefore has the local 
shape. We should note that this type of intrinsic NG 
(with local shape) also occurs in model of non-local in- 
flation where similar self- interactions are important [20| . 

In addition to the intrinsic NG, the non linear terms 
in Eq. (|3|) will also contribute to the bispectrum. In this 
case, Gaussian fluctuations in 8\ are transferred non- 
linearly to 8N, inducing non-Gaussianity. Because the 
non-linear terms in the SN expansion are local in space- 
time, they will give a contribution of the local form and 
the leading term is 



JNL ~ — _ ^3/2 



d X 3^4 



1 



7 2 k 2 



(12) 



where we have neglected higher order terms in the 5N 
expansion. Note that these two contributions have the 



same shape but opposite signs. The ratio of the intrinsic 
contribution to /jvl compared to the local contribution 



£ hit 
JNL 



f lo 



NL 



1 JlKxXX^ 2 (13) 

3 7)X H 2 1 j 



For the case where 7 is a constant (as in [7|), vanishes 
and the intrinsic contribution is the most important one. 
For the case where 7 is a linear function of Xj we have 
7,x ~ l/x an d the ratio can be simplified to 



P - r] x N eK 2 



(14) 



up to an order 1 coefficient which depends on the exact 
form of the potential. Given that one must have r\ x < 
N^ 1 in order for the damping (c. f. Eq. (O) to be small, 
this ratio is usually smaller than 1 although for values of 
rj x ~ 0.01, it can be of order 1. 

Finally there are higher orders (loop) contributions to 
the 2-pt and 3-pt that we have not included here. For 
example, from Eq ©, there will be a correction to the 
2-pt curvature function of the following form 



(a 2 ) 



loop 



H 



7,x7 ( <>Xk 



d 3 k' 
(2^ 



rSXk-k'SXk' 



(15) 

This is a loop factor and the integral is IR divergent with 
a log that must be regulated. A natural cutoff is to use a 
box of the size of the universe today 2l| . Then the loop 
factor is of order ~ j and using the intrinsic 3-pt 

function of Sx computed previously, we get a term of the 
order of 



%xl N e K 6 V, 



\n(kL) 
: 2vr 2 



(16) 



For a given model, it is possible for this term to dominate 
over the 7 2 term in Eq. although for the perturbation 
theory to remain under control, one should require that 
7 is small enough to insure that S(j> 2 is the dominant 
term in the 2-pt. In addition one can get an additional 
contribution to the 2-pt proportional to 



H 



7 2 x 



d 3 k'd 3 k" 

■— — ^— (Sxk-k>Sxk>Sxk~k>>Sxk>- 



(2tt) 6 



(17) 

This would give a contribution proportional to the square 
of the 2-pt function (of the order of e ^ M -2 7 2 x ^ 4 m kL) as 
well as a possible 2-loop contribution proportional to the 
4-pt function. If we assume that 7, xx is very small, we 
can stop here as all subsequent terms would be negligi- 
ble. A systematic approach to calculate all these loops 
contributions has been developed recently in 22j. Sim- 
ilar terms contribute to the 3-pt function as well [23j ]. 
A more thorough analysis of these additional contribu- 
tions is important and we leave a detailed calculation of 



4 



these terms for further work. For the example we discuss 
later, we parametrically estimated all the loop contribu- 
tions and chose a point in the parameter space where we 
can neglect the loops. Nevertheless, we need to point out 
that we also found region in the parameter space where 
these loops corrections must be included and can even 
dominate. 

Some comments are in order about the validity of 
the perturbative analysis. Assume for the moment that 
7 iX ~ such that only is non-zero. Then, in order 
to be able to compute the 2-pt correlation function of 
d\ using only the quadratic action like we did, the in- 



241 ] . In our context, it is 



teraction terms must be small 
enough to demand that V. xxx < H , or equivalently, 

/ fint \ V 3 

7>2(P|) 1 / 6 %^ , (18) 



where we used Eq. (jlip and we assumed e» ~ e/ and 
7 < 1. This tells us that in order to obtain a large Jnl 
we cannot take the coupling between the two sectors (x 
and 0) to be arbitrarily weak. Indeed, the smaller 7 cor- 
responds to a more strongly interacting scenario for a 
given /jvx- Note that this is a computational limit, it 
might be very interesting to have a strongly coupled hid- 
den sector. We leave this to future investigation. In the 
case where 7 x is non zero then one must ensures that 
the loop factors discussed above are small corrections to 
the 2-pt function. This can be achieved by taking 7 suf- 
ficiently small. 

A multi-field scenario. We illustrate this method of 
generating non-Gaussian curvature perturbations with a 
slightly modified version of a D-term inflationary sce- 
nario arising in IBMs This IBM scenario uses the 
inherent extra gauge symmetry of the hidden sector to 
reduce the amount of fine-tuning required for inflation, 
and can accommodate Standard Model embeddings with 
interesting phenomenology. We will find that this well- 
motivated scenario will exhibit the rich structure of non- 
Gaussian perturbations which we have discussed above. 

The key feature of this inflationary setup is that we 
have 4 different D-brane stacks with C/(l)i n f,i,2,3 gauge 
theories on their world volumes. There are several chiral 
multiplets which live at the topological intersections of 
any two branes with each other and transform in bifun- 
damental of the gauge groups on the two branes. The 
resulting D-term potential and superpotential are 

„2 



V D c = 
y inf 



>NG\ 



D 



9l 



(l<M 2 HS| 2 + 



1^ 



-0" 



-6) 2 



w 



q\S NG \ 2 ) 2 + - 
\S4>+4>- + \ngSng4>+4>ng 



(19) 



where • • • denotes additional fields which arise at the in- 
tersections of these 4 D-brane stacks with other branes 
which arc not relevant for us. V^ f is the D-term poten- 
tial of the brane which generates inflation, and q = ±1 is 
the charge of Sng under f7(l) 2 (in an IBM construction, 
this sign is determined by the orientation of topological 
intersection at which the multiplct lives). Note that we 
have made a particular choice of the sign of the charge 
for all other fields, but we have left the charge of Sng 
undetermined because we will find that it has interest- 
ing observational consequences. The superpotential will 
generate F-terms (setting = (f>N G = 0) 



V 1 



2 o2 1 2 



X'S 



+ x 



VG^NG' 



(20) 



Without any loss of generality, we assume £ > 0. As 
a result, 4>+ 1S the waterfall field. S is the inflaton; its 
vev gives <fi+ positive mass. Sng is the new field added 
compared to the set-up in [f| and this will be our field x- 



On the D-flat direction, V r f st 



and V& 



*R 2 . An 



inflaton potential is induced from the one-loop Coleman- 
Weinberg potential 



V ini = 
Vcw(x) = 



1 



9 



;V C w(x) 



9*e 
2 

(x 2 + l) 2 m(:r 2 + 1) - 2x 4 lnx 2 

+ {x 2 - l) 2 ln(a; 2 - 1) -4 In 2 , (21) 



where x — ^ and we have neglected x contributions 
to the 1-loop potential (in the regime we will consider, 
they are small). Using U{1) rotational invariance, we 
may define 



= ReS, x = ReS 



NG, 



T = Re< 



(22) 



and we can set all the imaginary parts to zero. The 

hidden sector field x nas a quartic coupling given by 

2 2 1 2 

Vhid = jx 4 = 92 2 93 X 4 - Finally the mass of the tachy- 
onic field depend on the quadratic sum of two fields 



m 2 T = -g 2 Z + >?4> 2 + (a?vg - 9.92k 2 



(23) 



Inflation ends when the inflaton reaches the critical value 

j,2 _ g 2 C 1 (ggl-Awc)x 2 

t^c A 2 A 2 



which gives 



dx 



ggj - X NG X 
A 2 6 r . 



(24) 



We denote the stochastic part of the field x by x, which 
is a random variable with a mean of zero. The quantum 
perturbations of x grows to be of order H at horizon exit. 
At this point, one can assume a semi-classical treatment 
where the field x is undergoing a random walk with steps 
of order H for each interval of time 1/H. The Fokker- 
Planck equation then gives the non-Gaussian probability 
distribution for x- We thus find x = o~H / y/v, where a 
typical range is —0.6 < a < 0.6 0]. Note that 7 depends 
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on this stochastic parameter and this is one of the main 
difference from [7J. 

Analysis. At this point the system has been reduced 
to D-term inflation, fn our analysis we will take the limit 



and 



In this limit 



and assume that 7 <C 1 as well as <f> 2 
- - ■ - 1 ~ — 'or equivalently e* 



B 1 
X 2 



h ~ ^j£- (or equivalently e* ~ e/), 
and the inflaton does not move much during inflation. 
We must go to this limit in order to realize significant 
non-Gaussianity because the perturbations to SN can 
arise both from fluctuations at the start of inflation and 
from the end. The initial fluctuations (in our model) arc 
Gaussian, so significant non-Gaussianity can only arise 
if the non-Gaussian perturbations at the end of inflation 
are of about the same scale. This also means that the 
spectral index will be close to one. In this regime, the 
2-pt. normalization is 



P c ~ 2.28 x 10" 



in 2 
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6 (21n(2)) 2 A 2 M p 6 



(25) 



The local and intrinsic contributions to /jvl are 
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A 2 £ 



(26) 



The sign of /at l is determined by the sign of qg 2 — X 2 NG . 
Note that this dependence arises from the way these cou- 
plings affect the mass of the waterfall field; the Yukawa 
coupling adds a positive mass contribution and delays 
the end of inflation, while the gauge coupling term adds 
a mass contribution whose sign depends on the charge q, 
and can either delay or hasten the end of inflation. In 
particular, if q = +1 and g\ > X NG , then qg\ — X 2 NG > 0. 
If this model arises from the type of IBM construction 
described in ||, then one would indeed expect g\ > Xfj G , 
because Yukawa couplings are exponentially suppressed 
(in Type IIA string theory, they arise from worldshect in- 
statons). On the other hand, if q = — 1, then we can set 
Xng = 0, because the Yukawa term XngSng'P+4'ng is 
not classically gauge-invariant (both <f>+ and Sng would 
have charge —1 under [7(1)2)- We then find that qg\ < 0. 
We thus directly see how the sign of the charge of the new 
scalar determines the sign of the various contributions to 
/nl in this model. 

There are a few consistency conditions one must sat- 
isfy; we take rj x < 0.01 so that the non-Gaussian pertur- 
bations in x do not decay significantly between horizon- 
crossing and the end of inflation. For perturbative con- 
trol we must satisfy Eq. (TT8"]) . As a toy example, we take 
C~2.9x 10- 6 Af 2 , A ~ 2 x 10- 4 , qg\ - X 2 NG ~ -0.17, 
v ~ 0.05, g 2 ~ 10~ 3 and \ ~ 0.12H/y/v. We then satisfy 
the above consistency conditions and generate 7 ~ —0.5 



Tint 
JNL 
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1.002 , 



45 



G/j~7x 10" 



(27) 



where we included the cosmic strings produced in such 
models (with Gfi — £/4). For this point in parameter- 
space, one can check that the loop corrections Eq. (fTS! 
[TBI scales as N e v 2 and H 2 /\ 2 respectively and for the 
parameters chosen above, they are small. But clearly 
for other points in parameter space one or both of these 
expansions parameters could be significant, in which case 
one would have to compute higher-order terms as well. 

As is typical in D-term inflation, this model predicts 
cosmic strings. When their contribution to the power 
spectrum is correctly accounted for, n,s ~ 1 can be com- 
pletely consistent with WMAP data [llj. One expects 
that larger values of £ can lead to smaller values of n s 
(though if the resulting cosmic string tension exceeds ob- 
servational bounds, the strings must be unstable). In 
this case, large values of /nl would require larger values 
°f 192 ~ ^ng- We can also find points in the parameter 
space with large f^L where Gfi is smaller (and satisfy the 
bound of [25|]) but in this case, the WMAP bound on n s 
while not excluding our model, does not favor it. Also, we 
find that the smaller the cosmic string tension the harder 
it becomes to satisfy all the consistency conditions and 
yet get large Jn l ■ We have assumed that the vev of the 
inflaton is sub-planckian. This implies that the energy 
density in gravity waves will be unobservably small. Sig- 
nificant gravity waves can be generated if the inflaton 
traverses super-Planckian scales, but in that limit the ef- 
fective field theory is no longer valid. Finally, we should 
mention that a different mechanism for generating NG 
from the tachyon was presented in (26|. 

Discussion. We have illustrated a simple model for 
generating primordial non-Gaussian curvature perturba- 
tions. This model can arise easily from intersecting brane 
models and allows a large /jv l with either sign which one 
expects to settle soon from the results of WMAP and 
Planck. The shape of the non-Gaussianity is local. One 
can fit the CMB data with large non-Gaussianity and a 
flat spectrum once the contribution from cosmic strings 
is included correctly. Without the cosmic strings, the 
current bound on n s implies this model is not favored by 
the data but not excluded. The fact that the tachyon 
mass depends on additional fields is a feature shared by 
many string theory models of inflation 271 ] . 

In this general and well-motivated effective field theory 
model, we have found a relatively complicated structure 
to the non-Gaussian contributions, with contributions 
from both intrinsic and non-linear terms, and possible 
loop-contributions as well. The non-linear contributions 
are usually dominant. We have found that the sign of 
Jn l depends largely on the hidden-sector charge of the 
extra field, and on the relative importance of the gauge 
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and Yukawa sectors. 

Another model of D-term inflation in string theory 
was recently reanalysed in 1 281] (and effects from entropy 
modes were discussed in [29J) and it will be interesting to 
study whether NG can also be generated in their model. 
An important difference in the D3/D7 model compared 
to the IBM model we presented is the fact that in the 
former, the Yukawa coupling and gauge coupling are the 
same. Therefore to achieve the n s ~ 1 regime, they need 
a very small gauge coupling (of order 10 -4 ) while in our 
case it can be larger. Of course, the D3/D7 model con- 
tains a full string theory description with backreaction 
and moduli stabilization effects included. It would be in- 
teresting to see if NG of the type presented in this work 
still survives in a more concrete model from string theory. 
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